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Abstract 

^ I A two-segment Fermi-Pasta-Ulam lattices has been investigated by using nonequilibrium molec- 

■ ular dynamics. Here we present an anomalous negative differential thermal resistance (NDTR) 

O . that have not been reported in Frenkel-Kontorova and 04 lattices up to the present. The NDTR 

disappears at low temperature region. The region of NDTR shifts from the large to the small 
' temperature difference region as the system size increases. Anomalous dependence of NDTR on 

(N : 

'nI" ■ the temperature can be explained as the negative effect induced by the nonlinear coupling. The 

0\ ■ 

! explanation can also cover the phenomenon of NDTR in momentum-nonconserved lattices. 
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I. INTRODUCTIONS 

Negative differential thermal resistance (NDTR), which may implies an intrinsic physical 
mechanism for some thermal devices, has attracted lots of attentions in recent years. NDTR 
is first found in asymmetric Frenkel-Kontorova (FK) lattices [1], and then in pure anhar- 
monic systems j^. Analytical and numerical studies have been carried to reveal the original 
property of NDTR [llQ. Up to now NDTR has been found in many low-dimensional sys- 
;ems, such as the anharmonic lattice with gradient mass two FK lattice with weak link 
l|, |5|, the pure 04 and FK lattice , the double-stranded systems jol, and so on. It has 
been reported that NDTR depends on the nonlinear external potential, the finite size and 
the temperature of the system. 

The main explanations about the phenomenon of NDTR include several points: 1) the 
mismatch of the phonon spectra of the two particles [l,^]; 2) the competition between the 
temperature difference, which acts as an external field, and the ternperature-dependent ther- 
mal boundary conductance I2I; 3) the effective phonon-band shifts [t]; 4) a ballistic transport 
that induce the competition between the molecular occupation factor and the temperature 
difference jsj; 5) the phonon-lattice scattering becomes so significant that NDTR occurs 



All these explanations suggest a common view that the NDTR cannot occur when 
the system size increases, or in the system without nonlinear external potential, or in the 
system as the temperature is increasing. In other words, NDTR only appears at low temper- 
ature in a small system with nonlinear external potential. For the sake of comparison, here 
we call this kind of NDTR the normal NDTR. According to the normal NDTR, it is pro- 
posed that the momentum-conserved system without nonlinear external potential, such as a 
Fermi-Pasta-Ulam (FPU) lattice, behaves no NDTR phenomenon. Since the nano-materials 



such as the carbon nanotubes is a momentum-conserved system [10|, , this view about no 
NDTR in the momentum- conserved system once places the studies of the NDTR in carbon 
nanotubes in a difficult position. Actually, it is not the case. 

In this paper, a two-segment FPU lattice with a weak link is studied by molecular dynam- 
ics simulations. We will present that the NDTR can also occur in a momentum-conserved 
system. It is surprising that the NDTR in FPU lattice behaves an anomalous dependence on 
the temperature and the system size, which is different from the NDTR found in the system 
with nonlinear external potential. Here we introduce a competition between the linearity 
and the nonlinearity to interpret the origin of NDTR. Our explanation can cover some of 
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the phenomena of NDTR that have been reported in FK models. 

II. MOMENTUM-CONSERVING AND MOMENTUM- 

NONCONSERVING LATTICES 

Momentum-conserving lattice [l2| A chain of N coupled atoms, in which only nearest- 
neighbor interactions will be considered for simplicity. The first class of models we wish to 
consider are defined by an Hamiltonian of the form {pi = rriixi) 
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Boundary conditions need also to be specified by defining xq and x^+i. Typical choices 
are periodic, fixed or free boundaries. As only internal forces, which depend on relative 
positions, are present, the total momentum is conserved and thus a zero mode exist. The 
important examples are the well-known Lennard- Jones potential and Fermi-Past a-Ulam 
(FPU) potential. 



Momentum-nonconserving lattice 12| At the simplest level of modelization, this can be 
described by adding an external, on-site potential to Eq.(l). For instance, neglecting the 
transverse motion leads to one- dimensional models of the form: 
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The substrate potential U{xi) breaks the invariance Xi — )■ + const: of Eq.(2) and 
the total momentum is no longer a constant of the motion. An important example is the 
well-known Frenkel-Kontorova (FK) potential. 

III. MODELS AND SIMULATIONS 

The nonlinear lattices that we use in this letter consist of two segments, left segment (L) 
and right segment (R). Each segment is a FPU lattice. Segment L and R are coupled via a 
spring of constant kint- The total Hamiltonian of the model is 



H = Hl + Hpi + Hint-, 



(3) 



and the Hamiltonian of each segment can be written as 
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with XM,i and XM,i denote the displacement from equihbrium position and the conjugate 
momentum of the z*'* particle in segment M, where M stands for L or R. The parameters /cm 
and Pm are the harmonic and anharmonic spring constant of the FPU lattice, respectively. 
We couple the last particle of segment L and R via a harmonic spring. Thus, Hint = 
%^ {xl,n — xr,nT ■ For the sake of simplicity, we set the mass of particle m = 1 and the 
Boltzmann constant ks = 1- 

In our simulations we use fixed boundary condition and the chain is connected to two 
heat baths at temperature Tl and Tr. It is reported that the Nose- Hoover thermostat may 
induce an inauthentic NDTR in the system without nonlinear external potential jo], lisl . 
To avoid the adverse effect of the Nose-Hoover thermostat we use the Langevin heat baths 
and integrate the equations of motion by using the Verlet frog-jumping algorithm HQ. 
The local temperature is defined as Tj = (pf). The local heat flux is defined as ji = 
kAiiPii^i — + f3M{Pi{xi — Xi^iY), and the total heat flux is J = Nj. The simulations 

are performed long enough to allow the system to reach a steady state in which the local heat 
flux is constant along the chain. For the sake of comparison, we define a heat current ratio, 
Jn = J / Jmax, in which Jmax is the maximum heat current under a fixed temperature Tr of 
the right heat bath. The transport coefficient is an important quantity for characterizing the 



transport mode of a thermal transport process [16|, The thermal conductance evaluated 
as cr = Nj/AT represents an effective transport coefficient that includes both boundary 
and bulk resistances 12|]. Here we fix the temperature of the hot heat baths and change the 
temperature of the cold heat baths. The temperature difference increases with the decrease 
of the temperature of cold heat baths. 

IV. RESULTS AND DISCUSSIONS 

Figure 1 shows the temperature dependence of NDTR in the two-segment FPU lattice. 
In the lattices with nonlinear external potential, which also called momentum-nonconserved 
systems, lots of studies have shown that the NDTR disappears with the temperature in- 
creasing of the systems j^, 5-7|. However, here it is presented that the NDTR occurs with 
the temperature increasing of the hot heat baths. Interestingly, the results for momentum- 
conserved systems are contrary to those for momentum-nonconserved systems. 

This anomalous NDTR can be understood from the competition between linear and 
nonlinear interactions. Firstly we would like to discuss the role of linear and nonlinear in- 
teraction on the thermal conductivity. It is well known that positive effects of the linear 
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FIG. 1: Heat current as a function of the temperature difference for various fixed high temperature 
Ti=0.05, 0.1, 0.35, 0.60, 1.0. The hnear and nonhnear couphngs are respectively k = 1.0 and 
/3 = 0.5, the weak hnk couphng A;,„j=0.01, and the system size is 64. 

interaction on the thermal conductivity of the system have been confirmed by many studies, 



which include not only numerical but also analytical results j^, ll2|, |18| . It can be described 
as a function K'X"", where K and A are respectively thermal conductivity and the linear 
coupling, the exponent a > 0. Correspondingly, it is also shown that the nonlinear external 

I . The relation- 
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potential preforms a negative influence on the thermal conductivity 
ship is written as K~V^, in which V is the strength of the nonlinear coupling or the external 
potential, the exponent /3 < 0. The NDTR depends exactly on the competition between 
these negative and the positive effects. When the negative effects surpass the positive ones, 
the NDTR occurs. Otherwise, the DNTR disappears. 

For the FPU lattices, figure2 shows that the heat current decreases with the nonlinear 
coupling constant, which implies a nonlinear inhibition on the thermal conductivity. In 
the absence of nonlinear coupling, i.e., in the case of a harmonic lattice, as shown in Fig. 
2, when the nonlinear coupling /3 equals zero, the heat current increases linearly with the 
temperature difference and no NDTR occurs. When /3 increases to a value larger than 0.3, 
due to the enhancement of the nonlinear effect the NDTR appears at the region of the large 
temperature difference. 

In order to give more details for this explanation, we plot figureS to show the dependence 
of the linear and nonlinear forces on the displacement. As mentioned above, in FPU lat- 
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FIG. 2: Heat current as a function of temperature difference for various nonlinear couplings. The 
temperature = 0.35, and the remain parameters are the same as those for Fig.l. 

tice the NDTR disappears for low temperature and occurs for high temperature. For the 
system with low temperature, the particles oscillates near the equilibrium position and the 
displacement between two nearest neighbor particle is small, i.e., AxfQ. As shown in Fig. 
3(a), w hen Axj~0, the force from linear coupling {—kAxi) is larger than that from nonlinear 
coupling (— /3Axf). Therefore, the linear coupling determinate the thermal transport and 
no NDTR occurs. As the temperature of the system increases, Axi\ goes to a value larger 
than Ax^l and the force from — /3Axf is larger than —kAxi. The nonlinear coupling as the 
negative effect determinate the thermal transport and then NDTR occurs. However, in the 
FK lattices an opposite temperature dependence of NDTR is observed. The Hamiltonian 



of FK model is Hpj, = Eti [S^ + f - XM,^f + ^(1 - cos(^a;,)) 

Pi denote the displacement from equilibrium position and the conjugate momentum of the 
jth poj-f^ioiQ Tj^g parameters /cq and V are the harmonic spring constant and the strength 
of the external potential, respectively. Here we can also give a explanation from the rela- 
tionship between the force and the displacement. As shown in Fig. 3(b), for the system 
with low temperature, the force from linear coupling {—koAxi, ko = 0.8) is smaller than 
that from nonlinear coupling (— V^sinxj, V = 0.6) and the nonlinear coupling determinate 
the thermal transport and then NDTR occurs. On the contrary, for the FK lattices with 
high temperature, the superior linear coupling induces the disappearance of the NDTR. If 
one set ko = 3.0, the linear force is always larger than the nonlinear force, then the NDTR 



„2 



, where Xj and 
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disappears, which have been confirmed in Refs. 
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FIG. 3: Force depends on the displacement of the particle, (a) linear and nonlinear coupling force 
in FPU model, (b) linear coupling and nonlinear on-site force in FK model. 

Figure 4 shows how the coupling constant hint affects the temperature dependence of 
the heat flux J/Jmax- The region of NDTR diminishes for increasing until it vanishes 
at a critical coupling constant. The weak link hint is a linear coupling, as also depicted 
in Fig. 3 (a), when hint is small, the nonlinear effect can easily get above the linear one at 
Axi~0 in the interface of the two FPU lattices and then NDTR occurs. Instead, the NDTR 
vanishes for strong linear link hint- 
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FIG. 4: Heat current as a function of temperature difference for different weak link couplings ki^t- 
Here the remain parameters are k=1.0, /3=0.5, Ti;,=0.35 and N=64, respectively. 
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In addition, the system size dependence of NDTR in FPU lattices is also different from 
that have been found in FK and 04 lattices. In FK and 04 lattices {2], Q , with the increasing 
of the system size, the NDTR region shifts from the small to the large temperature difference 
region and vanishes in the end. On the contrary, as shown in Fig. 5, in FPU lattices the 
NDTR disappears for the small system. When the system size increases, the region of NDTR 
(the regime with a frame) shifts from the large to the small temperature difference region. 




FIG. 5: NDTR depends on the system size. Here the remain parameters are k=1.0, /3=0.5, T/,=0.35 
and kint=0-05, respectively. 

V. COMPARISONS 

Here we would like to compare our explanations about NDTR with those have been 
presented. As mentioned above, the NDTR have been explained as the phonon-lattice 
scattering, the crossover of the transport mode, the competition between the temperature 
difference and the boundary conductance, the mismatch of phonon spectra and so on. As 
have been reported, the mismatch of phonon spectra and the effective phonon-band shifts 
can give a good interpretation for the NDTR in two-segment lattices, but fail to interpret the 
NDTR in pure single lattice . Moreover, the competition between the temperature 

difference and the temperature-dependent thermal boundary conductance can only explain 
the small systems with strong boundary effect j^, ^, but will be invalid for large systems 
mentioned here. Finally, the phonon-lattice scattering and the crossover of the transport 
mode approach the main reason of NDTR at low temperature j^, 5, 6| , but not an appropriate 
origin of NDTR at high temperature. Actually, all these explanations are based upon a 



8 



negative effect, which maybe imply a fundamental physical mechanism of NDTR. And here 
this negative effect comes from the nonlinearity that includes both the nonlinear coupling 
and the nonlinear external potential. 

A comparison of NDTR in the two cases is also presented from three aspects: 1) the 
temperature dependence of NDTR. In Ref. 5|], a temperature dependence of NDTR is 
reported. When the temperature of hot heat baths Tr is small, there exists NDTR. As Tr 
increases, the NDTR disappears. However, here it is presented that the NDTR occurs with 
the temperature increasing of the hot heat baths; 2) the system size dependence of NDTR. 
In the previous works (FK and phi_4 model), as shown in Refs. {2!, Sill, NDTR becomes 
weak and disappears with the increasing of the system size. The most important property 
is the NDTR region shifts from the small to the large temperature difference region and 
vanishes in the end. In the current works, on the contrary, as shown in Fig. 5, in FPU 
lattices the NDTR disappears for the small system. When the system size increases, the 
region of NDTR (the regime with a frame) shifts from the large to the small temperature 
difference region. 3) the nonlinear potential. As shown in Eqs. (1) and (2), the Hamiltonian 
of momentum-conserving lattice is much different from that of momentum-nonconserving 
lattice. Many reports have shown that the momentum-conserving lattice has an anomalous 
thermal conductivity 19(|. Then it is naturally to regard the NDTR of momentum-conserving 
lattice as an 'anomalous' one. Additionally, the NDTR in this paper performs very different 
or even opposite behavior from the NDTR that have been reported in the previous papers, 
then we have to call it 'anomalous NDTR'. 

VI. CONCLUSIONS AND APPLICATIONS 

In summary, we have reported an anomalous negative differential thermal resistance in 
two-segment FPU chains through molecular dynamics simulations. The NDTR occurs in the 
system when the system size as well as the temperature increases. This anomalous NDTR 
effect can be understood from the negative effect induced by the nonlinear force, which 
can give a more fundamental explanation about NDTR. Since FPU lattice is a momentum- 
conserved system, our results provide a new view about the NDTR in the system without 
external potential, which is regarded as a necessary condition for NDTR up to the present. 
Furthermore, our results have also suggested that the NDTR can be achieved in two-segment 
nanoscale materials with a weak linear link, which maybe a exciting information for fabri- 
cating a nanoscale thermal device. 



9 



The experimental motivation of the model refers to the carbon nanotubes and their ex- 
tended models. The Lennard- Jones potential or Brenner-Tersoff potential, which describes 



momentum-conserving lattice, is usually used to study carbon nanotubes [lO|, lUl- A ad- 
vanced structure of the nanotubes is more similar to our model. For example, one can 
connect two carbon nanotubes with silicon nanotubes. Here silicon nanotubes are weak 
connections and carbon nanotubes are momentum-conserving lattices. 
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